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ABSTRACT 


Multiplicative  Rayleigh  fading  is  a  frequent  problem  in  wireless  commimications.  If  the 
channel  is  relatively  benign  and  fading  is  not  severe,  one  may  obtain  higher  bit  rates  for 
an  equivalent  bandwidth  by  using  M-ary  Quadrature  Amplitude  modulation  (MQAM). 
A  variation,  used  to  combat  channel  fading  while  still  retaining  MQAM,  is  differential 
MQAM  (i.e.,  DQAM).  The  term  differential  refers  to  the  phase  that  is  coded  exactly  as  in 
DPSK;  however,  the  amplitude  is  still  subject  to  distortion  by  the  fading  channel.  In  this 
paper,  we  propose  a  technique  called  quotient  coding,  which  is  designed  to  remove  channel 
effects  from  the  symbol  amplitude  as  well  as  its  phase.  In  particular,  we  shall  apply 
it  to  MQAM,  resulting  in  modulation  which  we  term  Quotient  Quadrature  Amplitude 
Modulation  (QQAM).  In  contrast  to  DQAM,  QQAM  is  just  as  effective  at  suppressing  the 
effects  of  channel  fading  with  respect  to  the  entire  symbol  as  DPSK  is  for  the  phase  alone. 
Furthermore,  it  can  be  applied  to  arbitrary  MQAM  configurations. 
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1  INTRODUCTION 


Differential  coding  presents  a  standard  technique  for  combating  channel  instability- 
in  wireless  communications  [l]-[7].  Moreover,  differentially  encoded  signals  do  not  require 
synchronization  for  suppressed  carriers,  thereby  avoiding  additional  noise  from  devices  such 
as  phase- locked  loops.  This  approach  is  particularly  effective  when  the  basic  modulation  is 
M-ary  phase-shift  keying  (MPSK).  in  which  case  it  is  labeled  DPSK  (differential  MPSK). 
On  the  other  hand,  if  the  channel  is  relatively  benign  and  fading  is  not  an  issue,  one  may 
obtain  higher  bit  rates  (as  well  as  an  improved  signal-to-noise  ratio  (SNR)  performance 
over  MPSK)  for  an  equivalent  bandwidth  by  using  M-ary  quadrature  amplitude  modulation 
(MQAM).  A  variation,  particularly  used  to  combat  multiplicative  Rayleigh  fading  while 
still  retaining  MQAM,  is  differential  MQAM  (i.e.,  DQAM).  In  the  case  of  DQAM,  the 
differential  encoding  removes,  or  at  least  mitigates,  the  effects  of  channel  fading  upon  the 
phase  component  of  the  MQAM  signal.  However,  the  amplitude  is  still  vulnerable  to  the 
channel,  and  one  must  resort  to  approximate  techniques  that  estimate  the  magnitude  of 
fading  (e.g.,  automatic  gain  control  (AGC))  in  order  to  extract  the  symbol  amplitude  [3].^ 
Such  methods  fail  if  the  fade  rate  becomes  too  large.  That  is,  they  suffer  from  a  trade-off 
between  the  quality  of  the  channel  estimation  and  the  speed  at  which  it  is  tracked  [7]. 
In  this  paper,  we  propose  a  technique  called  quotient  coding,  which  is  designed  to  remove 
channel  effects  from  the  symbol  amplitude  as  weU  as  its  phase.  In  particular,  we  shall 
apply  it  to  MQAM,  resulting  in  a  modulation  that  we  refer  to  as  Quotient  Quadrature 
Amplitude  Modulation  (QQAM).  Unlike  DQAM,  QQAM  is  just  as  effective  at  suppressing 
the  effects  of  channel  fading  with  respect  to  the  entire  symbol  as  DPSK  is  for  the  phase 
alone.  Furthermore,  it  can  be  applied  to  arbitrary  MQAM  configurations. 

Our  discussion  takes  place  at  baseband.  In  fact,  since  pulse  shape  is  irrelevant  to  the 
discussion,  the  description  shall  be  entirely  in  terms  of  complex  symbols.  To  develop  our 
notation  and  set  the  context,  we  begin  with  a  brief  outline  of  MPSK,  MQAM,  and  their 
differential  counterparts.  MPSK  modulation  may  be  summarized  by  [5]: 

5  =  ^e2xj(m-i)/M  ^  (la) 

where  A  is  real, 

(l)m  =  27r;(m-l)/M;  m  =  1,  . . . ,  M  ,  (16) 

and  M  is  the  total  number  of  symbols  to  be  transmitted. 

We  shall  assume  that  the  channel  in  question  can  be  modeled  by  multiplicative  fading 
(typically  Rayleigh  [4])  and  additive  noise.  In  other  words,  the  received  signal  takes  the 
form  a{t)s{t)+n(t),  where  a{t)  is  a  complex  number  with  appropriate  time  correlation,  and 
n(t)  is  noise.  Thus,  for  example,  to  demodulate  MPSK  either  one  must  estimate  the  phase 
of  o,  which  is  difficult  when  the  fading  is  rapid,  or  one  must  modify  the  modulation  scheme. 

^Star  16QAM  differentially  encodes  1  bit  of  amplitude,  but  the  technique  does  not  generalize  to  2  or 
more  amplitude  bits  [l]-[3],  [6],  [7]. 
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A  standard  technique  is  to  use  DPSK  (differential  MPSK),  which  may  be  implemented  as 
follows:  Let  d(t)  be  the  transmitted  DPSK  signal,  and  define  its  amplitude  Ad(t)  and 
phase  9(t)  by 

dit)  =  (2) 

Note  that  Ad{t)  is  real.  Assume  that  d(0)  is  known,  for  example  that  it  equals  1.  Then, 
define  d{t)  as  a  function  of  its  previous  value,  d{t  —  1),  and  of  the  current  symbol,  s{t),  by 

d(t)  =  Ad(t)  =  ^d{t-l).  (3) 

Let  dr  be  the  received  signal,  and  a  “bar”  indicate  complex  conjugation.  The  phase  of  the 
symbol  s(t)  is  recovered  by 


phase(  s(t)  )  =  phase(  dr{t)  dr{t  —  1)  ),  (4) 

and  is  essentially  independent  of  the  fading  channel.  That  is,  it  is  independent  to  the 
extent  that  the  phase  shift  due  to  the  channel,  phase(^  ^(0  —  1)  }■>  is  approximately 

zero. 

MQAM  is  a  combination  of  amplitude  and  phase  modulation.  That  is, 

s  =  Ar,  =  Ane^^p  n=l,  ...,iV;  p=l,  ...,P,  (5) 

where  the  total  number  of  symbols,  M,  is  given  by  M  =  NP.  DQAM  is  implemented  by 
differentially  encoding  the  phase: 


and  recovering  the  symbol  by 


Sr(t)  -  <*-(')  Jill  ■  (66) 

But  only  the  phase  of  Sr(t)  is  impervious  to  the  channel.  The  received  amplitude,  |5r(t)|  = 
|Q;(t)s(t)|  requires  channel  estimation  in  order  to  extract  the  amplitude  of  s(t).  This  may 
be  done  by  using  averaging  techniques  such  as  AGC,  but  there  is  a  performance  trade-off 
inasmuch  as  one  must  average  (a)  enough  to  smooth  out  the  dependence  on  A„,  but  (b) 
over  a  short  enough  time  period  to  track  the  channel. 
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2  QUOTIENT  CODING 


The  basic  modulation  paradigm  for  quotient  coding  is  quite  simple.  One  modulates 
via  the  relationship 

q{t)  =  s{t)q{t-l),  (7a) 

and  demodulates  by 

Sr(t)  =  -  1).  (76) 

(The  tildes  appear  in  equations  (7)  to  distinguish  q  from  our  final  result  in  equations  (10) 
and  (13)  below.)  It  is  clear  that  the  multiplicative  channel  cancels  out  both  with  respect 
to  phase  and  amplitude.  That  is,  in  the  absence  of  noise, 

qr(t)  /  qr(t  -  1)  =  (a(t)q(t))  /  (a(t  -  l)q(i  -  1))  «  q(t)  /  q(t  -  1).  (8) 

There  is,  however,  a  practical  difficulty  associated  with  equation  (7a).  If  |s(t)|  ^  1,  the 
transmitted  power  will  experience  large  excursions  caused  by  the  variation  in  the  amplitude 
of  the  transmitted  symbol.  In  particular,  \q{t)\  =  |  nt.<f  l9(0)l-  On  the  average, 

these  excursions  will  average  out.  But  there  are  still  likely  to  be  long  sequences  producing 
very  small  or  very  large  transmitted  powers.  The  former  is  sensitive  to  noise,  the  latter 
physically  unrealizable.  (The  2DASK  differential  amplitude  scheme,  which  recognizes  a 
“1”  by  a  change  in  amplitude  and  a  “0”  otherwise,  avoids  this  problem  if  it  is  implemented 
so  that  the  transmitted  symbol  takes  on  exactly  two  amplitudes;  but  this  is  only  possible 
for  N  =  2.  When  used  in  conjunction  with  8PSK,  such  a  scheme  is  generally  referred  to 
as  Star  DQAM  [1],  [2],  [6],  [7].)  This  dynamic  range  problem  is  alleviated  by  a  simple,  but 
critical,  modification  which  we  proceed  to  describe. 

Let  us  return,  momentarily,  to  the  MQAM  modulation  of  equation  (5).  Without  loss 
of  generality,  we  order  the  symbols  by  size, 

Ai  <  A2...  <  An,  (9a) 

and  normalize  the  amplitudes  so  that 


Ai  =  1.  (96) 

Let  Pmax  be  the  maximum  transmission  power  at  baseband.  The  modified  quotient  coder 

is  defined  by 

r  q(t  -  1)  s(t)  |3(t  -  1)  s(t)|  < 

9(0  =  <  _  (10) 

(  q(t  —  1)  /  s(t)  other. 

A  reasonable  initialization  value  is  q{-l)  =  y/P^/2.0.  Now,  let  Pr(t)  be  the  instanta¬ 
neous  received  power, 

Pr(.t)  =  l9,-(i)l".  (11) 
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Note  that  in  the  absence  of  noise,  Pr(t)  /  Pr{t  -  1)  is  independent  of  the  fading  a{t),  and 
that  it  equals  |5(t)p  or  1  /  |5(<)p  according  to  equation  (10).  Since  Ai  =  1  >  1  /A^ ,  these 
two  ratios  correspond  to  two  populations:  that  for  which  y/Pr{t)  /  Pr(t  -  1)  is  greater  than 
or  equal  to  Ai  and  that  for  which  it  is  less  than  or  equal  to  I/A2.  (For  the  case  |s|  =  Ai, 
the  ratios  of  the  powers  are  equal;  i.e.,  |s(t)p  =  l/|s(t)p  =  1.) 

Define 

,  ^  0.6(4,  +  i) 

=  0.5(1  +  j-y  (12) 

The  demodulation  scheme  becomes 

(  qr{t)  /  qr{t  -  1)  for  Pr{t)  >  7]^  Pr(t  -  1) 

=  \  _  _  ,  (13) 

(  qr{t  -  1)  /  qr{t)  for  Pr{t)  <  Tj'^  Pr(t  -  1) 

which  may  also  be  written 


I  qrit)  qr{t  -  1)  /  Prit  -  1)  for  Pr{t)  >  Pr(t  -  1) 

Srit)  =  <  _  .  (14) 

(  qr{t)  qr(t  -  1)  /  Prit)  for  Prit)  <  Prit  -  1) 

It  should  be  remarked  that  the  actual  constellation  transmitted,  equation  (10),  is  a  non- 
linearly  distorted  version  of  the  associated  MQAM  constellation  equation  (5).  However, 
it  does  closely  resemble  MQAM  locally,  over  the  transmission  of  a  few  symbols  (cf.  Ap¬ 
pendix).  In  fact,  the  detector  behavior  is  roughly  equivalent  to  an  MQAM  constellation 
under  a  mixture  of  of  SNRs.  In  addition,  one  must  take  into  account  the  probability  of  an 
error  in  the  binary  decision  based  on  the  threshold,  (cf.,  equation  (13)). 
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3  AN  EXAMPLE 


To  illustrate,  we  present  an  example  with  two  amplitude  bits  and  three  phase  bits;  in 
other  words,  32QQAM  with  iV  =  4  and  P  =  8.  The  constellation  on  which  this  example 
is  based  is  pictured  in  figure  1.  Its  rings  are  equally  spaced  with  amplitudes,  Ai,  given 
by  the  formula  Ai  =  1  +  0.4  (i  —  1),  and  are  individually  Gray  coded.  The  factor  0.4  was 
obtained  by  running  several  simulations  and  picking  the  best  value.  However,  we  have  not 
made  any  serious  attempt  to  optimize  the  constellation.  There  is  no  reason  to  assume  that 
equally  spaced  rings  or  that  annuli  are  the  best  geometry.  For  that  matter,  what  is  “best” 
depends  on  the  application. 


Figure  1.  Underlying  QQAM  (and  also  the  DQAM)  constellation  used  in  the  simulation  example. 
There  are  3  bits  assigned  to  phase  and  2  bits  to  amplitude. 


The  DQAM  modulation  used  for  comparison  was  also  based  on  the  constellation  of 
figure  1.  To  be  fair,  several  runs  were  made  to  determine  whether  the  ring  locations 
should  be  adjusted,  but  it  was  found  that  the  DQAM  behavior  was  relatively  insensitive 
to  the  ring  spacing  and  that  0.4  was  close  to  optimal  under  equal  spacing.  An  AGC  with  an 
exponential  window  was  used  to  handle  fading.  More  precisely,  the  channel  was  tracked  via 
Pav{ti)  =  (1  ~  Pavi^i-i)  +  A  |<ir(t)P,  where  Pav  represents  the  average  power  estimate. 
A  value  of  A  =  1/40  =  0.025  was  chosen.  It  corresponds  to  a  window  of  40  symbols.  This 
value  was  originally  determined  by  optimizing  performance  for  the  fade  rate  under  study 
(0.0018  normalized  Doppler^  ),  but  was  found  to  be  close  to  optimal  even  for  very  slow 
fading  (fade  rates  of  10”®  and,  perhaps,  lower).  Roughly  speaking,  it  is  about  as  short  a 
window  as  one  can  use  while  still  averaging  sufficiently  over  the  four  possible  amplitudes. 

Figure  2  compares  32QQAM  with  32DQAM  in  a  fading  environment  with  no  additive 
noise  (Ej/iVo  =  oo).  The  advantage  of  QQAM  is  rather  remarkable.  Except  at  extremely 
low  fade  rates  (cf.,  figure  4),  the  AGC-dependent  DQAM  suffers  much  more  from  the 
fading  than  QQAM  does.  This  huge  performance  gap  is  somewhat  more  subdued  in  the 


■^The  normalized  maximum  Doppler  frequency  is  defined  as  max-doppler-frequency/symbol-rate  = 

fd/ f symbol  ■ 
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presence  of  noise.  Figure  3  contains  simulation  results  as  a  function  of  the  signal-to-noise 
per  bit  Eb/No  for  a  fading  channel  where  the  maximum  normalized  Doppler  frequency  is 
fixed  at  0.0018.  At  high  SNR,  where  fading  is  the  dominant  effect,  QQAM  does  much 
better  than  DQAM.  More  precisely,  there  is  a  crossover  point  at  about  12  dB.  Past  this 
point,  32DQAM  saturates,  and  32QQAM  gains  dramatically  in  performance. 

The  other  extreme  from  figure  2,  no  fading  at  all,  is  illustrated  in  figure  4.  For 
completeness,  we  have  also  included  performance  curves  for  32PSK,  32DPSK,  and  32QAM. 
Since  there  is  no  fading,  this  case  is  effectively  to  the  left  of  the  crossover  point,  and 
32DQAM  significantly  outperforms  32QQAM.  The  gap  is  rather  large,  but  this  is  not 
surprising  inasmuch  as  the  32DQAM  of  figure  4  was  implemented  by  assuming  perfect 
knowledge  of  the  channel  amplitude. 


Normalized  doppler  freq  (fd  /  fsymb) 


Figure  2.  Comparison  of  32QQAM  and  32DQAM  with  3  bits  of  phase  and  2  bits  of  amplitude. 
The  channel  model  consists  of  pure  multiplicative  Rayleigh  fading  [4]  (i.e.,  no  additive  noise),  and 
the  x-axis  is  expressed  as  the  normalized  maximum  Doppler  frequency  (i.e.,  fd/ fsV'^hol). 
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Figure  3,  Comparison  of  32QQAM  and  32DQAM  as  a  function  of  EhfNo*  The  channel  model 
consists  of  multiplicative  Rayleigh  fading  with  a  normalized  Doppler  frequency  of  0.0018  plus 
additivie  white  gaussian  noise  (AWGN). 


Figure  4.  Performance  of  32QQAM  and  several  other  modulation  schemes  in  the  presence  of 
pure  AWGN.  No  AGO  was  used  for  the  32DQAM  (i.e.,  it  was  frozen  at  the  true  power).  From 
left  to  right,  the  curves  are  32QAM,  32DQAM,  32PSK,  32QQAM,  and  32DPSK. 


4  CONCLUSION 


We  have  described  a  very  general  technique  for  removing  the  deleterious  effects  of  mul¬ 
tiplicative  fading  from  most  modulation  schemes.  Termed  quotient  coding,  it  determines 
the  value  to  be  transmitted  by  forming  the  quotient  (or  the  product,  when  dictated  by 
the  current  state)  of  the  previous  value  with  the  current  symbol.  Quotient  coded  MPSK 
is  equivalent  to  DPSK.  The  application  of  quotient  coding  to  an  annular  16QAM  con¬ 
stellation  produces  a  variant  of  Star  DQAM  [6].  Quotient  coding  of  an  arbitrary  MQAM 
constellation  results  in  a  modulation  form  that  we  have  named  QQAM.  The  actual  constel¬ 
lation  transmitted  is  a  nonlinear,  distorted  version  of  the  associated  MQAM  constellation. 
Nevertheless,  it  resembles  MQAM  locally,  i.e.,  over  the  transmission  of  a  few  symbols.  In 
fact,  the  QQAM  detector  behavior  is  roughly  equivalent  to  an  MQAM  constellation  under 
a  mixture  of  of  SNRs. 

The  advantage  of  QQAM  is  that  it  removes  the  effects  of  channel  fading  from  both 
the  amplitude  and  phase.  Differentially  coded  modulations  such  as  DQAM  still  leave  the 
amplitude  bits  vulnerable  to  fading,  a  problem  that  must  be  dealt  with  by  equalization 
techniques  such  as  AGC.^  As  an  example,  we  have  examined  the  relative  performances  of 
32QQAM  and  32DQAM.  We  found  that  at  low  Eb/No,  where  the  effects  of  additive  noise 
dominate  the  fading,  DQAM  performs  better  than  QQAM,  but,  at  high  Ei/No,  DQAM  is 
fading-limited  and  falls  behind.  For  the  case  studied  (0.0018  normalized  maximum  Doppler 
frequency)  there  was  a  crossover  point  at  about  12  dB.  Past  this  point,  32DQAM  saturates, 
and,  as  a  result,  QQAM  dramatically  (by  several  orders  of  magnitude)  outperforms  DQAM. 

The  intention  of  this  paper  was  to  introduce  a  modulation  technique  which,  as  in¬ 
dicated  by  simulation,  has  distinct  advantages  in  fading  environments.  Many  questions 
remain.  A  full  analytical  characterization  of  QQAM  performance  would  be  interesting, 
but  appears  difficult.  We  have  not  broached  the  behavior  of  non- annular  geometries  or  the 
optimization  of  performance  relative  to  geometric  parameters.  A  mathematical  treatment 
of  the  distribution  of  the  magnitudes  of  the  transmitted  values  could  also  provide  some 
insights. 


o  ^ 

An  exception,  Star  DQAM,  applies  differential  coding  to  the  amplitude,  but  the  method  is  limited 
to  a  single  amplitude  bit. 
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APPENDIX:  Constellations  and  the  Computation  of  BEP 


A  complete  theoretical  treatment  of  the  bit  error  probability  (BEP)  for  QQAM  is 
well  beyond  the  scope  of  this  paper/study.  However,  it  is  enlightening  to  at  least  derive 
some  general  formulae  and  outline  the  relationship  of  QQAM  to  the  underlying  MQAM 
Constellation.  When  noise  is  present,  the  demodulation  takes  the  form  (cf.,  equation  (13)), 


qi  +  ni 
?2  +  n2’ 


(A.1) 


where  qi  and  52  represent  q{t)  or  its  conjugate  at  two  successive  times,  and  ni  and  n2,  the 
corresponding  noise.  Assuming  large  SNR,  we  may  make  the  approximation. 


gi  +  rii  _  gi/g2  +  ni/q2 

g2  +  n.2  1  +  ^2  /52 

?2  g2  g2  g2  J 

«  il  _  !^ 

g2  \  gl  g2 


(A2) 


Clearly,  the  factor  91/92  has  no  effect  on  the  signal-to-noise  ratio.  In  fact,  letting  cr^ 
be  the  noise  variance,  and  a  and  b  be  the  magnitudes  of  91  and  92,  respectively,  (A.2) 
implies 


7(0,  b)  =  effective  SNR/ symbol 

1 

(aVa2)  +  (cr2/62) 

1  a2fe2 


<t2  a2  +  62 

Es 

=  -T  6^. 
cr2 

(A.3) 

where 

Es  =  P(a2)  =  E(62), 

(A.4) 

and 

=  E,  a^  +  IP- 

(A5) 

Note  that  /i  is  independent  of  the  scaling  of  a  and  6;  i.e.,  of  any  scaling  of  the  signal.  The 
bit  error  probability,  Pb,  can  be  computed  by  first  conditioning  on  a  and  b.  That  is. 


P qqam  —  P(a,  b). 

a^b 


(A.6) 


10 


One  may  consider  approximating  Ps^a^b)  by  (a)  computing  the  probability  of  error 
under  the  effective  SNR  which  is  ^(a,  b){Eb/No)  and  (b)  by  taking  into  the  the  probability 
of  an  erroneous  decision  with  respect  to  rj  of  equation  (13).  However,  summation  over  a 
and  b  is  much  more  formidable  than  one  would  expect.  With  very  few  exceptions  (e.g.,  N 
=  2),  the  values  of  q{t)  generated  via  equation  (10)  form  an  infinite  Markov  chain.  More 
precisely,  suppose  that  N  =  2,  A2  =  2,  and  Pmax  =  4.  Then,  the  possible  transitions  from 
q(t-l)  to  q(t)  are  that  1  may  go  either  to  1  or  2,  and  2  may  go  to  2  or  1.  This  is  illustrated 
in  figure  5a.  However,  for  example,  if  N  =  3,  then  the  Ai  are  1,  2,  3,  and  y/P^  =  3,  a 
potential  sequence  of  transitions  starting  with  q  =  2.0  is  2,  2/3,  4/3,  8/3,  8/9,  16/9, 16/27, 
32/27,  . . .  corresponding  to  transmitted  symbols  3,  2,  2,  3,  2,  3,  2,  ...  (cf.,  figure  5b). 
Ultimately,  q{t)  takes  on  an  infinite  number  of  values. 


2 


(b) 

Figure  5.  Illustration  of  propagation  and,  hence,  of  the  distribution  of  the  values  q{t).  The 
small  numbers  are  the  symbol  magnitudes  responsible  for  the  transition.  The  nodes  are  the  values 
of  q(t).  In  figure  5(a),  N  =  2  and  MQAM  amplitudes,  Ai  =  1  and  A2  =  2.  In  figure  5(b),  N  =  Z 
and  MQAM  amplitudes,  Ai  =  1,  A2  =  2,  and  A3  =  3. 

A  consequence  of  this  is  that  the  actual  received  QQAM  constellation  differs  con- 
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siderably  from  the  underlying  MQAM  configuration  For  example,  figure  6  is  the  received 
constellation  (i.e.,  the  eye  diagram)  for  32QQAM  with  no  noise  and  no  channel  fading 
over  5000  symbols.  The  constellation  fills-in  with  respect  to  amplitude  (i.e.,  the  spokes). 
In  practice,  this  poses  no  difficulty  for  the  detector  which  really  relies  on  a  ratio  of  am¬ 
plitudes  which,  int  turn,  have  the  distribution  of  the  underlying  MQAM  constellation.  In 
fact  we  see  from  equations  (A.2)  to  (A. 6)  that,  locally,  the  detector  behaves  very  simi¬ 
larly  to  MQAM  under  a  mixture  of  SNRs.  On  the  other  hand,  analytically  computing  its 
performance  is  a  nontrivial  problem. 


Irmginary  Pan 
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Figure  6.  Constellation  for  32QQAM  (no  noise,  no  fading,  5000  symbols  plotted). 
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